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Centrality evolution of pt and yt spectra from Au-Au collisions at ^/snn = 200 GeV 
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A two-component analysis of spectra to pt = 12 GeV/c for identified pions and protons from 200 
GeV Au-Au collisions is presented. The method is similar to an analysis of the rich dependence of 
Pt spectra from p-p collisions at 200 GeV, but applied to Au-Au centrality dependence. The soft- 
component reference is a Levy distribution on transverse mass mj. The hard-component reference 
is a Gaussian on yt with exponential {pt power-law) tail. Deviations of data from the reference are 
described by hard-component ratio vaa which generalizes nuclear modification factor Raa- The 
analysis suggests that centrality evolution of pion and proton spectra is dominated by changes in 
parton fragmentation. The structure of vaa suggests that parton energy loss produces a negative 
boost Ai/t of a large fraction (but not all) of the minimum-bias fragment distribution, and that 
lower-energy partons suffer relatively less energy loss, possibly due to color screening. The analysis 
also suggests that the anomalous p/iv ratio may be due to differences in the parton energy- loss 
process experienced by the two hadron species. This analysis provides no evidence for radial flow. 

PACS numbers: 13.66.Bc, 13.87.-a, 13.87.Fh, 12.38.Qk, 25.40.Ep, 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

Pt spectrum analysis of relativistic heavy ion colli- 
sions typically invokes the assumption that a thermally- 
equilibrated flowing bulk medium is produced. The hy- 
drodynamic (hydro) model plays a central role in the de- 
scription [U, [2, Q . Hard-scattered partons are proposed 
as probes of the medium A thermalized system is 
said to expand and cool until chemical and kinetic decou- 
pling. Pt spectra for identified hadrons are used to derive 
thermodynamic and chemical properties of the medium. 
Blast-wave fits and mass dependence are used to infer ra- 
dial flow Temperatures and chemical potentials are 
inferred from statistical-model spectrum fits and particle 
ratios Q ■ Unexpected spectrum properties motivate sup- 
plemental recombination or coalescence models [8, 9, 10]. 
From the combined spectrum information evidence for a 
QCD phase transition and formation of a quark-gluon 
plasma is sought. 

Observationally, pt spectra for identified hadrons from 
heavy ion collisions at RHIC extending to pt = 12 GeV/c 
as in this analysis may vary over 10 orders of magnitude, 
whereas relative data accuracy may be 5-10% over much 
of the spectrum. It is thus difficult to represent visu- 
ally the full information content of carefully measured 
spectra. According to conventional assumptions pt is 
separated into a thermal/hydro region pt < 2 GeV/c, 
a coalescence/recombination (ReCo) region 2 < pt < 6 
GeV/c and a perturbative QCD (pQCD) region pt > & 
GeV/c. Each interval receives a separate analysis strat- 
egy and theoretical treatment. One can question whether 
all physical information in the spectra is being used to 
test theory. In this paper I show that differential spec- 
trum analysis is essential to exploit the data fully. 

Analysis of the multiplicity dependence of pt spectra 
from NSD p-p collisions revealed that the two-component 
model of pt spectra provides an accurate and com- 
plete description of p-p spectrum data for unidentified 



hadrons jll| . The soft-component shape is independent 
of multiplicity. The hard component shape is indepen- 
dent of multiplicity, but its amplitude relative to the soft 
component is proportional to multiplicity. The p-p hard 
component is approximately a Gaussian on transverse 
rapidity yt, and its form is related both to parton mo- 
mentum spectra and to parton fragmentation functions 
in e+-e~ collisions p^ . 

The hard component of p-p spectra can be interpreted 
as a minimum-bias parton fragment distribution (FD) — 
conditional fragmentation functions folded with the unbi- 
ased parton spectrum. In that interpretation a substan- 
tial fraction of the full hadron spectrum down to pt 0.3 
GeV/c is part of the parton fragment distribution, with 
implications for A-A pt spectra. Parton fragmentation 
then competes with hydro, thermalization and recombi- 
nation to account for the bulk of A-A particle production. 
To sort out competing mechanisms we must understand 
A-A Pt spectra to the statistical limits of the data, with 
as little physical model dependence as possible. 

In this paper I apply the two-component spectrum 
model to pion and proton spectra from five centr ality 
classes of Au-Au coUisions at ^snn — 200 GeV T3|. 
Spectra are compared to a model function, extended from 
the two-component analysis of p-p collisions, which de- 
scribes all structure in the data. 

The paper is organized as follows. I introduce the data 
used and discuss some conventional aspects of spectrum 
analysis. I describe the two-component spectrum model 
as applied to p-p and Au-Au collisions. I discuss differ- 
ential spectrum analysis, including conventional nuclear 
modiflcation factor Raa and its generalization in this 
analysis to hard-component ratio rAA- I also introduce 
a simple parton energy-loss model including attenuation 
and negative boost of the minimum-bias fragment distri- 
bution. 

I apply the two-component model to pion and proton 
spectra, obtaining soft- and hard-component model func- 
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tions for each species and five centralities. I then obtain 
hard-component ratios vaa for the data which generaUze 
Raa and Rcp results. I model the data vaa trends with 
simple functional forms which suggest that parton energy 
loss produces a negative boost of part of the minimum- 
bias fragment distribution (the reference hard compo- 
nent). I compare those results to theory and observe 
that the data suggests the presence of color screening 
for lower-energy partons, explaining the large abundance 
of minijets observed with correlation analysis in central 
Au-Au collisions. 

Finally, I demonstrate that the anomalous p/ir ratio 
may be an aspect of parton energy loss differently man- 
ifested for pions and protons, compare integrals of the 
spectrum models to the Kharzeev and Nardi (K-N) two- 
component model [14], and obtain centrality trends of 
(pt) which are compared to published values. 



II. CONVENTIONAL pt SPECTRUM ANALYSIS 

The invariant single-particle density on 3D momentum 
space for identified hadrons at ry = is 



p{xt) = l/27r 1/xt (fin/dxt di] 



(1) 



averaged over 2t: azimuth and one unit of pseudorapidity, 
with transverse variable xt = pt, "nit or yt for momentum, 
mass or rapidity respectively. The Jacobian for transfor- 
mation from Pt to yt is pt rrit /yt , that for pt rrit is 
unity. Although the data analyzed here are for identi- 
fied hadrons I use pseudorapidity to provide convenient 
comparisons with unidentified-hadron spectra. 

Ratio techniques have been introduced to study "jet 
quenching" [isj . The conventional ratio measures (nu- 
clear modification factors) are 
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the latter applied to more peripheral (P) collisions rel- 
ative to more central (C) collisions when p-p reference 
spectra are unavailable for comparison. The ratio limits 
for Pt small are l/v or vpjvc-, where mean participant 
pathlength v = 2nunary / nparUdpant [Hi- Some features 
of the spectra at smaller pt are suppressed by those defi- 
nitions. One purpose of this analysis is to restore access 
to those features. 



A. 200 GeV Au-Au PID pt spectra 

The particle-identified (PID) data used in this analy- 
sis are large-statistics identified pion and proton pt spec- 
tra extending to pt ~12 GeV/c for five Au-Au centrali- 
ties . Spectra were obtained from 15M 0-12% central 
Au-Au collision events and 14M minimum-bias events at 



^/sWn = 200 GeV. The centralities are 0-12%, 10-20%, 
20-40%, 40-60% and 60-80%, corresponding to mean par- 
ticipant pathlength v = 5.50, 4.87, 3.93, 2.83 and 1.93 re- 
spectively [l^ . Spectra are in the form of invariant par- 
ticle density p{pt) — p{xt) in Eq. ([T]) with xt ^ Pt- 



B. Assumptions and interpretations 

The analysis in (isj represents a conventional pt spec- 
trum analysis. By hypothesis spectra are divided into 
three pt intervals: 1) Particle production in < 2 GeV 
is said to be dominated by soft QCD processes and de- 
scribed by hydro models assuming local thermalization 
and collective flow [ll,i,[3]. 2) Inp* G [2, 6] GeV/c differ- 
ing meson and baryon suppression is interpreted in terms 
of hadronization via "coalescence of constituent quarks 
from a collective partonic system" [l3|. 3) For 

Pt > 6 GeV/c particle production is assumed to reflect 
pQCD parton scattering, parton energy loss and frag- 
mentation. 

Two Pt intervals are relevant to this analysis. 

1) The intermediate-pt region: The stated issue is 
relative suppression of protons and pions. Protons are 
said to be "less suppressed" than pions, but significant 
suppression is still observed, in contrast to d-Au colli- 
sions where significant enhancement (Cronin effect) is 
observed [13]. Differences between proton and pion sup- 
pression (anomalous p/tt ratio) motivate the quark coa- 
lescence or recombination model, predicated on the prior 
existence of a thermalized collective partonic medium 
(QGP) |,i[i|. 

2) The large-pt region: Hadron yields arc strongly sup- 
pressed in this Pt interval, and pQCD theory attempts to 
predict the suppression in terms of parton energy loss. 
A major issue is whether suppression of different hadron 
species can reveal color or flavor dependence of parton 
energy loss [H, [ll|. Pion and proton suppression are 
similar, and the p/tt ratio returns to that for d-Au, sug- 
gesting that collectivity and coalescence are not impor- 
tant in this interval. 



III. TWO-COMPONENT SPECTRUM MODEL 

The two-component spectrum model is the exten- 
sion of a model of hadron production in nuclear colli- 
sions [3j[20'] . To describe pt spectra hadrons are assumed 
to emerge from two processes, soft and hard, distinguish- 
able by parametric variations (e.g., centrality v or rich)- 
Each process has a characteristic functional form on pt 
or yt which is approximately independent of collision pa- 
rameters, but the relative amplitude of the two processes 
varies smoothly with a control parameter, providing the 
basis for separation. In the present paper the analysis is 
treated as a model-independent Taylor expansion of den- 
sity p{yt',v) on V, with no physical model assumptions. 
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The terms of the series can then be identified with the 
conventional physics-based two-component model. 



A. Transverse rapidity yt 

pt spectra provide information about parton fragmen- 
tation, but the fragmentation process is more visually ac- 
cessible on a logarithmic momentum variable p^ . The 
dominant hadron (pion) is ultrarelativistic at 1 GeV/c, 
and a relativistic kinematic variable is a reasonable alter- 
native. Transverse rapidity in a longitudinally comoving 
frame near midrapidity yz — is defined by 



yt = ln([mt +pt]/mo), 



(3) 



with mo a hadron mass and = + uiq. p{yt) is de- 
fined by Xt yt in Eq. ([J) ; transformation of data den- 
sities p{pt) to p{yt) employs the Jacobian below Eq. ([T]). 
For pions (corresponding to all yt values in this pa- 
per) Pt — 0.16, 0.5, 1, 2, 4, 6, 10 GeV/c transforms to 
yt = 1, 2, 2.66, 3.33, 4, 4.5, 5 respectively. 

The initial version of this analysis defined separate yts 
for pions and protons according to their masses (toq) 
to accommodate a radial flow component according to 
the Cooper- Frye formalism ^21,] , with a common positive 
boost Ayt (blue shift) of nearly-thermal spectra on trans- 
verse rapidity. However, there is no evidence for hydro 
structure in the spectra, and the hard component for pro- 
tons was found to be approximately the same on pt as 
that for pions. Thus, yt was subsequently defined strictly 
as a logarithmic transformation from pt using the pion 
mass for both pions and protons. The transformation 
increases visual access to spectrum structure but has no 
other physical significance. For thermal production from 
a boosted source the proper yt for each hadron would 
be required, and spectra should be considered in general 
with both yt definitions in mind. 



B. p-p two-component model 



Event multiplicity rich serves in effect as a hard trigger. 
P = rich ~ 2.5 is the pt-integrated ID NSD hadron den- 
sity on Tj averaged over 27r azimuth, and r/ e [—0.5,0.5] 
at 200 GeV. Because the mean minijet multiplicity and 
dn/drj for NSD p-p collisions are both ^ 2.5, x is ap- 
proximately the probability to observe one hard parton 
scatter in -q G [—0.5,0.5] in one NSD p-p collision. 

Soft-component model function Sq is a unit-normal 
Levy distribution 



1 



mt - ntp 



(5) 



an apparently minor modification (m,t vs pt) of the 
"power-law" distribution A/(l -I- pt/po)" @, but repre- 
senting "soft" rather than "hard" processes (longitudinal 
rather than transverse fragmentation) [ll[. The density 
on vfit is defined with the proper pion or proton mass, 
but the transformation to S^iyt) is then made with the 
same pion mt and yt used for both protons and pions. 

Hard-component model function _ffo in [Tl| is a unit- 
normal Gaussian on yt 



Ho{yt;yt,(Tyt) = ^^exp 



-{yt - ytf 



(6) 



As with 5*0 the form is independent of control parameter 
rich- In this Au-Au analysis the common yt with pion 
mass is used for protons and pions, but the Gaussian 
parameters are found to be somewhat different for the 
two hadron species. 

For this analysis po = d^n/d'qd4) — 2.5/27r is the 
average pt-integrated 2D hadron density on (ry, (p) at 
mid-rapidity. The pt-integrated soft and hard compo- 
nents are p^ — po/{l + x) and ph — pox/{l + x), with 
X = 0.012 for NSD p-p coUisions at 200 GeV. I then de- 
fine Sppiyt) = PsSoiyt), Hpp{yt) = phHo{yt), and the 
multiplicity-dependent per-particle 3D density is 

Ps/Ps ■ Ppp{yt;nch) = Sppiyt) + rich/nch ■ Hpp{yt). (7) 

That expression for p-p collisions is formally equivalent 
to Eq. ^ for A-A coUisions. 



In the two-component spectrum model applied to p- 
p collisions at 200 GeV the soft-component density (~ 
longitudinal nucleon fragmentation) is simply propor- 
tional to the integrated event multiplicity Uch within 
T] £ [—0.5,0.5] [rich dn/drj = p, the ID density on 
77), and the hard-component density transverse par- 
ton fragmentation) is proportional to the multiplicity 
squared [llj . yt spectra from p-p collisions are accurately 
described by Ppp{yt) = Ps{nch)Soiyt) + ph{nch)Hn{yt), 
with Phi Ps = X ^ 0.005 Tie?! and ph + Ps = p {ps and ph 
replace and Uh from [ll|). Sq is a Levy distribution 
on mt and is a Gaussian on yt, each normalized to 
unity. In the notation of the present analysis 

PppiVux) = p{So{mt[yt]) + X HQ{yt)}/{l + x). (4) 



C. Ho and the power- law tail 

The spectra used for the two-component analysis 
in [11] extended only to pt ^ 6 GeV/c. The principal 
result was isolation of the hard component of the p-p 
Pt spectrum as a Gaussian on yt describing parton frag- 
ment distributions [l2| . In the present analysis the spec- 
tra extend to pt = 12 GeV/c, and the form of Hq must 
be modified. The QCD power-law trend, implemented 
previously by a monolithic "power-law" function [23j . 
is implemented in this analysis within modified hard- 
component model function Hq. The method to generate 
Hq as a Gaussian on yt with exponential tail consistent 
with pQCD power law p^"'' is described in App. [X] 
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The pQCD exponent Uh associated with the underly- 
ing parton pt spectrum should be clearly distinguished 
from the model parameter rig associated with the soft- 
component Levy distribution. In Fig. 2 (right panel) 
of [LI] the variation of "power-law" fitting parameter n 
with p-p event multiplicity ranges from Ug E [14, 20] for 
small multiplicities to Uh € [7, 8] for large multiplicities: 
the hard component increasingly dominates the spec- 
trum with increasing multiplicity (c/. Fig. 10, left panel 
of [HI). The conventional "power-law" function mod- 
els two physical mechanisms with one functional form. 
Its exponent n is oversubscribed, representing one pro- 
cess at smaller pt or rich and another process at larger 
Pt or rich- In the present analysis we observe Uh ~ 7.5, 
consistent with QCD expectations. Because of the Ja- 
cobian between pt and yt the exponential parameter in 
Ho{yt',yt, o-y^jUy^) is Uy^ ~ 5.5, as described in App. [Xl 



D. A- A two-component model 

The two-component model of hadron production (spec- 
trum integrals) in A-A collisions assumes that the soft 
component is proportional to the participant pair num- 
ber (linear superposition of N-N collisions), and the hard 
component is proportional to the number of N-N binary 
collisions (parton scattering) [l3|- Any deviations from 
the model are specific to A-A collisions and may reveal 
properties of an A-A medium. In terms of mean partic- 
ipant path length v = 2nhin/npart the pt-integrated 
Au-Au hadron density on rj is 



2 dn 

'^part dfj 



= Ps+ V Ph, 

= po [1 + 2^KN - 1)], 



(8) 



with skn = p/i/po and xkn ^ 0.08 for 130 GeV and 
0.10 for 200 GeV Au-Au colhsions [11]. From the p-p 
two-component analysis we expect xkn x — 0.012, 
which seems to conflict with [IJ]. I resolve the apparent 
contradiction in Sec. I XII 

In this analysis I extend the two-component particle- 
production model by analogy with Eq. ^ to describe 
Au-Au single-particle spectra on pt and yt- As with the 
p-p analysis the model is essentially the first terms of a 
Taylor expansion on A-A centrality parameter v 



(9) 



SNNiut) + vrAAiuuv) HNNiVt)- 



Ideally, the A-A reference functions would be p-p limit- 
ing cases Spp and Hpp. But there are small differences 
depending on what hadron species or combinations are 
used to define Xpp and X^n, and possible isospin effects. 
A-A limiting cases are therefore specified in terms of N-N 
collisions. 

With those exceptions the shape of Snn ~ Spp ap- 
pears not to change significantly with A-A centrality. 



For peripheral collisions we expect Haa Hj^n- For 
more central collisions Haa is strongly modified, reveal- 
ing parton energy-loss systematics which are absorbed 
into factor taa, the ratio of Haa to i/jvTV- 



IV. DIFFERENTIAL SPECTRUM ANALYSIS 

Conventional methods of differential spectrum analy- 
sis emphasizing the hard component at large pt are com- 
pared with alternative methods which reveal structure 
over the entire pt/yt acceptance. 



A. Conventional ratio measures 

Conventional spectrum ratios are expressed in two- 
component form as 



Raa 



Snn + V Haa 
Spp ~t~ Hpp 



(10) 



Rcp — — 



vp Sc + vc He 
vc Sp + vp Hp 



Effects of the medium on hard processes are measured 
at large pt by the deviation of Raa from unity according 
to ^] . The alternative form 
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(11) 



should go asymptotically to 1 at small pt or yt and to v 
at large pt or yt. I compare data in both forms to their 
two-component references. Since soft components 5* are 
approximately equal the limits of i?* for pt small, where 
dominates, are l/u or vp/i/c- However, important 
hard structure may persist at small pt which is strongly 
suppressed by the i?* definitions. 



B. Isolating the hard component 

To extract the hard component from data I define the 
difference 



Haa 



1 



-pAA — Snn 



(12) 



As noted, I adopt Snn ^ Spp from [11[ but optimize the 
model parameters for Au-Au data and the specific hadron 
species, so that 2/npart ■ PAAiut) Snn &s v ^ 0. The 
difference reveals hard component Haa for all pt or yt 
within the acceptance — all parton fragments, not just the 
pQCD region of the spectrum. In turn, I expect Haa — *■ 
Hnn as v ^ I. The limiting cases provide consistency 
checks on the analysis. An analogous procedure with the 
limit hch — was applied in flT] to isolate Spp and Hpp. 
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C. Parton energy- loss model 

A simple two-parameter model of the effect of parton 
energy loss on the hard component is 

HAAivuv) = A{v)HNN[yt + ^Vt{v)l (13) 

with A < 1 representing uniform attenuation and Aj/f < 
representing a negative boost or red shift on yt (uniform 
fractional momentum reduction) . I then define logarith- 
mic difference 

A \og{HAA{yt ; -)} = log ( ^^(^ I (14) 

as a generalization of ratio Raa which in contrast re- 
veals details of fragmentation systematics over all pt or 
yt- Differential analysis involves comparison of data to 
reference Eq. (|13p . Access to fragment distribution ratios 
over the entire yt acceptance, as illustrated in Sec. IVIIi 
should improve understanding of parton energy loss and 
fragmentation. 



D. Plotting formats 

Because of their great numerical range and strong cur- 
vature Pt spectra extending to large pt values don't pro- 
vide sufficient visual resolution to reveal all significant 
physical information. Dot sets displaced by powers of 
ten strongly suppress differential structure. More infor- 
mation is revealed by improved plotting strategies and 
differential comparisons with a two-component reference. 

Transverse rapidity yt has three advantages over pt for 
a differential analysis: 1) Parton fragment distributions 
have a simple, nearly-symmetric form on yt ^12]. 2) The 
Jacobian for pt — > yt removes a factor pf at larger yt , re- 
ducing the vertical range in this analysis by nearly three 
orders of magnitude and improving sensitivity to differ- 
ential structure. 3) QCD power law at larger pt 
transforms to a linear trend in a semilog plot on yt. 

Spectra are plotted as curves (line segments connect- 
ing data points) with common factor n^artl'^ removed to 
form per-participant-pair densities. For participant scal- 
ing such curves would coincide for all centralities. Varia- 
tions unique to heavy ion collisions appear as differences 
between the curves. Only plots with Haa isolated are 
sufficiently differential to reveal detailed structure. In 
those plots each centrality is identified by a specific line 
style. In other plots only the general relation of data 
(solid curves) to the two-component reference is shown. 

The dotted model curves in each plotting space provide 
a reference grid describing the two-component model. 
Plotting scales are coordinated between equivalent plots 
for pions and protons to aid comparison. This analy- 
sis reveals simple systematic deviations of data from the 
model. Correspondence of spectrum variations across the 
full yt acceptance are made apparent for the first time. 



V. TWO-COMPONENT ANALYSIS - PIONS 

Pt spectra for identified charged pions from five central- 
ities of 200 GeV Au-Au and hadron spectra from NSD 
p-p collisions are compared in several plotting formats, 
and the hard components are isolated. Au-Au spectra 
were renormalized from the published data as described 
in App. [Bl The spectra from ix'^ and tt" were combined, 
since ± ratios show no structure on pt . 

A. Pion spectra and soft component 

In Fig. [1] Pt spectra for pions from five Au-Au cen- 
tralities are plotted in the conventional format (five solid 
curves) [l^. The solid dots are the h^h = 1 spectrum 
from 3M NSD p-p collisions [11]. N-N hard component 
i^NN (dash-dot curve) is defined in Sec. lVBi including a 
power-law tail as described in App. |^ Spp and H^p (the 
latter not shown) from [ll| are redefined based on this 
analysis and combined to form the solid curve passing 
through the p-p points. 




; hydro ReCo pQCI) 



p^ (GeV/c) 

FIG. 1: Pion pt spectra for five Au-Au centralities (solid 
curves). The dotted curves are corresponding two-component 
model functions, bounded by dashed curves for N-N and 
6 = Au-Au collisions. The dash-dot curve is the N-N hard- 
component reference function. The dotted curves labeled Spp 
and Sjvjv are soft-component models for p-p and Au-Au colli- 
sions. The points are data from NSD p-p collisions (the solid 
curve through points is the two-component p-p model). A 
power-law trend p^^^ i^ shown at the right side of the plot. 

Soft component Satat (labeled dotted curve) has the 
same Levy form as for p-p collisions. The Au-Au 
pion parameters are = 20.2 ± 0.1, T = 0.1445 ± 
0.001 GeV, and = 12.0 ± 0.2. Exponent n, = 14.5 
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with As = 21.1 is used for the p-p hadron spectrum in 
this analysis. The p-p exponent obtained in [ll| was 
Us = 12.8. However, that analysis described hadrons us- 
ing a Gaussian hard component with no power-law tail. 
In the present analysis with power-law tail an improved 
description of the p-p data is achieved, and the "third 
component" described in [ll| is eliminated. It is clear 
from Fig. [T]that below pt = Q GeV/c {yt = 4.5) the true 
power-law trend is not easily detected. 

The dotted curves intersecting the data show the two- 
component reference for five values of v corresponding to 
the data centralities. The dashed curves bounding the 
dotted curves represent limiting cases v — 1, 6 for N-N 
and 6 = Au-Au collisions respectively. The power-law 
trend Pt~^'^ is illustrated. Suppression of the more-central 
spectra at large pt relative to the reference is qualitatively 
apparent. Otherwise, the detailed structure of individual 
spectra is better studied with more differential formats. 
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FIG. 2: Pion yt spectra for five Au-Au centralities (solid 
curves). The features are the same as for Fig. [1] The power- 
law trend e~^'^^' is indicated at the right side of the plot. 

In Fig. [5] the same spectra are plotted on yt- The 
spectra are nearly flat at small yt, aiding extrapola- 
tion and normalization. Hard-component reference Hj^n 
(dash-dot curve) is a Gaussian with exponential tail; the 
smooth transition to the power-law trend is apparent. 
The point of transition to the power law (in this case 
at yt — 3.75) depends on the centroid and width of the 
Gaussian and the power-law exponent. The plot on yt 
better displays the main part of the fragment distribu- 
tion centered near yt = 2.66 (pt ^ I GeV/c). 

In Fig.[3]i?yiA from data (differing line styles) is plotted 
on pt (conventional format) and yt for comparison. The 
thin solid curves are the two-component references for 
five centralities. The main inference from Raa on pt 




Pt yt 



FIG. 3: Left panel: Conventional nuclear modification fac- 
tor Raa for pious and five Au-Au centralities (thicker curves 
with changing line styles). The thin reference curves are ob- 
tained from the two-component model. Right panel: The 
same curves transformed to yt- The Raa limit for small yt is 
l/v, with 1/ the mean participant path length 



(left panel) is "jet quenching" — reduction of data from 
unity in the pQCD region > 6 GeV/c. Adopting unity 
as a reference [13] overestimates the magnitude of the 
reduction, since there is still a substantial soft-component 
contribution in that interval. The correct reference for 
each centrality is the corresponding thin model curve. 

Turning to Raa plotted on yt (right panel), the data 
for 60-80% (upper, bold solid curve) follow the reference 
within statistics, while more-central data fall below the 
reference on the right, but exceed the reference on the left, 
below Pt ^ 2 GeV/c. In a conventional description that 
part of the spectrum is also described as "suppressed." 
The format in the left panel, with unity taken as a(n in- 
correct) reference, gives a false impression of parton en- 
ergy loss and subsequent fragmentation. The Raa defini- 
tion suppresses hard-component structure in the small-pt 
region. The large-pt region and "jet quenching" are em- 
phasized at the expense of possible new physics at smaller 
yt ■ The entire fragment distribution should be compared 
to the correct two-component reference. 



B. Pion hard component 

In Fig. m the v Haa for pions from five Au-Au cen- 
tralities (bold curves of different styles) are plotted with 
the hard component from NSD p-p collisions [ll| (solid 
points). A common soft component Snn has been sub- 
tracted from the Au-Au data as in Eq. (fT^ . N-N hard- 
component reference H^n (lower, dash-dot model curve) 
is adopted from JdJ with essentially the same Gaussian 
parameters {Ah = 0.333 ± 0.005, yt = 2.66 ± 0.02, dy, = 
0.45 ± 0.005). However, the hard component for this 
analysis also includes a power-law tail with exponential 
constant Uy^. — 5.5 ± 0.3 (cf. App. The upper, 

dashed model curve is limiting case 6 Hj^ n corresponding 
to 6 = Au-Au. The dotted reference curves between de- 
scribe v HpfN for the five data centralities, the reference 
system. 
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FIG. 4: The hard component of pion yt spectra in the form 
vHaa (thicker curves with changing line style) compared to 
two-component reference z/ ffjvjv (dotted curves). The dashed 
reference curves are limiting cases for v = 1, 6. 



The agreement between Hnn ^ Hpp and the solid p-p 
points reveals the quality of the two-component repre- 
sentation for elementary collisions, consistent within sta- 
tistical errors. ffjvAr is also statistically consistent with 
the 60-80% Au-Au data for ~ 2. In fact, the value of 
Uy^ was determined by p-p and peripheral Au-Au data 
beyond yt = 4: {pt ~ i GeV/c). Incorporating the power- 
law tail improves the two-component description of p-p 
spectra. 

The remarkable new feature of this figure is the large 
excesses compared to the N-N reference at small yt for 
more central collisions and the apparent persistence of 
the fragment distribution there. Whereas there is a sub- 
stantial reduction of fragments at large yt there is strong 
enhancement at small yt , hinted in Fig. [3] (right panel) 
despite the l/v suppression. That pt interval is conven- 
tionally claimed for hydro and blast-wave fitting models. 
However, the trend with increasing centrality is a shift 
to smaller yt, not larger as one would expect for hydro 
expansion (positive yt boost). The inconsistency with 
hydro becomes more evident in the proton spectra. 



VI. TWO-COMPONENT ANALYSIS - 
PROTONS 

Identified pt spectra for protons from five centralities 
of 200 GeV Au-Au are shown in several plotting formats, 
and the hard components are isolated. The spectra were 
renormalized from the published data as described in 
App. E] The spectra for antiprotons are similar to those 



for protons, but there are two significant differences, de- 
scribed in Sec. lVIB] There are major differences between 
pion and proton spectra, visually apparent even for pt 
spectra. Differences below 6 GeV/c are conventionally 
attributed to a combination of hydro (radial flow) and 
quark recombination (anomalous p/tt ratio). I examine 
those interpretations in the context of the two-component 
model in Sec. |Xl 



A. Proton spectra and soft component 

In Fig.O proton spectra (solid curves) from ^Ts'l are 
plotted in the conventional format. The dotted curves are 
two-component reference curves for the five values of 
from data, bracketed by dashed curves for limiting cases 
v = 1, 6. Soft component Sjsin is a Levy distribution 
on proton mj. The proton parameters are = 3.87 ± 
0.05, T = 0.224 ± 0.005 GeV, = 17 ± 0.2, consistent 
with a previous analysis of proton spectra within pt G 
[0.3,3] GeV/c 




p, (GeV/c) 



FIG. 5: Proton pt spectra for five Au-Au centralities (solid 
curves). The general features are comparable to Fig. [1] 

As noted in Sec. IIII A( a substantial radial flow con- 
tribution was expected as part of the soft component, 
described by pion and proton Maxwell-Boltzmann dis- 
tributions on yt, with mass toq specific to each species 
but having a common positive rapidity boost Ayt (blast- 
wave model). However, the ly dependence for yt < 3.3 
{pt < 2 GeV/c) reveals that variations with centrality are 
consistent with a fixed soft component nearly Maxwell- 
Boltzmann in form and a hard component of fixed form 
approximately the same as that for pions. 

In Fig. [S] the same spectra are plotted on yt (solid 
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FIG. 6: Proton yt spectra for five Au-Au centrafities (soUd 
curves). The general features are comparable to Fig. [2] 



curves). Again, the spectra vary slowly near yt = 2, 
facilitating normalization and extrapolation to pt = 0. 
The soft component Snn is the labeled dotted curve. 
The other dotted curves show the two-component refer- 
ence for five values of v matching the data, bounded by 
dashed curves for v — I and 6. At lower right a power- 
law trend with Uy^ = 5.0 is sketched. Dash-dot curve 
if AT AT is the hard-component reference. 

In Fig. [7] (left panel) conventional ratio measure Raa 
for protons is plotted on yt . The apparent large suppres- 
sion near yt = 5 is similar to that for pions. The data 
in the region at intermediate pt ^ 2.5 GeV/c {yt ~ 3.5) 
are said to be "less suppressed" than those at larger pt, 
indicating that the reference value is assumed to be 1 
everywhere. The thin model curves are two-component 
references for the five centralities, and the actual refer- 
ence value for central collisions at ?/t ~ 3.5 is Raa ~ 0.3. 
Note the agreement within statistics between the 60-80% 
data (z/ ~ 2, top- most solid data curve) and the corre- 
sponding two-component reference (top-most thin refer- 
ence curve), similar to Fig. [31 confirming the soft and 
hard reference parameters. 

In Fig. [7] (right panel) v Raa goes asymptotically to h' 
at large yt- It is even more apparent in that plot that 
the peak near yt = 3.5 is a large excess relative to the 
correct reference, not a suppression. A possible origin of 
the proton peak, the source of the so-called "p/tt puz- 
zle," is discussed in Sec. |lXl Below yt — 2.66 the data 
agree well with the two-component reference, suggesting 
that parton energy loss may be significantly reduced for 
smaller parton energies. 




yt yt 



FIG. 7: Left panel: Conventional nuclear modification factor 
Raa for protons and five Au-Au centralities on yt (thicker 
curves with changing line styles). The thin reference curves 
are obtained from the two-component model. The Raa limit 
for small yt is Right panel: Alternative nuclear modifica- 
tion factor ly Raa compared to thin two-component reference 
curves. The limit for large yt is v. The most prominent fea- 
ture is a large excess near = 3.5 (pt ~ 2.5 GeV/c). 



B. Proton hard component 

In Fig. [8] I show vHaa for protons in five Au-Au 
centralities (bold curves with different line styles), and 
the hard-component reference Hj^n (dash-dot curve) in- 
ferred from analysis of the Au-Au spectra. There was no 
a priori model for the proton hard component. It was 
determined iteratively for this analysis just as for the p- 
p analysis. Spectrum variation with u was extrapolated 
to 1/ = to obtain the soft-component Levy parameters. 
The v variation of the remainder suggested the presence 
of a hard model function similar to that for pions and 
hadrons. 

The proton hard-component Gaussian model is de- 
fined by Ah = 0.475 ± 0.02, yt = 2.66 ± 0.05 and 
(jy^ = 0.33 ± 0.05, the width substantially smaller than 
the pion value 0.45. The power-law index Uy^ = 5.0 ±0.3 
appears to be significantly smaller (two sigma) for pro- 
tons than for pions (5.5). Because of the reduced Gaus- 
sian width the proton power-law trend starts earlier, at 
yt — 3.25. The hard-component reference with QCD 
power law thus plays an essential role from below 1 up to 
12 GeV/c, a much broader extent than the usual pQCD 
description applied above 6 GeV/c. 

The dotted curves represent the two-component refer- 
ence for five values of u. It is notable that the maximum 
values of the Hnn for protons and pions near yt — 2.66 
are nearly equal. Because the proton FD is narrower 
the proton fragment density becomes smaller at larger 
yt- The nearly equivalent shapes of the FDs may re- 
flect fragmentation to different hadron species from sim- 
ilar underlying parton spectra, terminating just above 1 
GeV/c. 

Spectra for antiprotons differ from those for protons 
by a change in the exponent Uy^ from 5.0 for protons to 
5.5 for antiprotons (the same as for pions) and a multi- 
plicative factor 0.8 relative to proton spectra (the factor 
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FIG. 8: The hard component of proton yt spectra in the form 
vHaa (thicker curves with changing line style) compared to 
two-component reference _ffjvjv (dotted curves). The general 
features are comparable to Fig. |31 



consistent with |25|). 

Au-Au proton data for all centralities are consistent 
with the two-component model below ~ 3 {pt ~ 1.5 
GeV/c), providing even stronger evidence (than pions) 
that radial flow does not play a substantial role in spec- 
tra. The major new feature is the large excess at yt 3.5 
for more-central collisions, which is discussed in Sec. IIXI 



plotting format is compatible with Eq. (|15p. The line 
rAA = 1 describes the two-component reference for all 
centralities (all dotted v H^n curves in previous plots). 
Deviations from 1 represent all residuals from the two- 
component reference, p-p data are consistent with the 
model for all rich- The structure in this figure is then 
unique to heavy ion collisions. The figure can be com- 
pared with Raa trends in Fig. [3] (right panel) . 
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FIG. 9: Hard-component ratios for pions and five Au-Au cen- 
tralities (thicker curves with changing line styles) relative to 
the N-N hard-component reference. The connected dots are 
data from NSD p-p collisions. The dash-dot lines represent a 
simple parton energy-loss scenario (cf. Sec. IVIlH . 



VII. HARD-COMPONENT RATIOS 

Hard-component references and data are compared di- 
rectly in the form of ratio taa = Haa/ H^jf which gener- 
alizes Raa as a measure of parton energy loss. Whereas 
Raa is called the "nuclear modification factor," taa can 
be called simply the hard-component ratio. Raa sup- 
pression at large pt is ambiguous because either an atten- 
uation factor A < 1 or a negative rapidity shift Ayt or 
some combination within the "pQCD" yt interval could 
result in Raa ^ 1 for > 6 GeV/c. In this and the 
next section I show that the ambiguity can be reduced 
by examining the structure of taa over the entire pt or 
yt acceptance. 

A. Pion ratios 

Fig.[9]shows Haa/Hnn vs yt, with Haa for pions from 
five Au-Au centrality classes (curves of several line types) 
and hadrons from NSD p-p collisions (solid dots). H^n 
is the model function defined in Sec. IV Bl The log-ratio 



Ratio data for peripheral collisions (60-80%) are also 
consistent with the two-component reference. For the 
other four centrality classes the trend above yt ~ 4.5 is 
consistent with Raa measurements in that interval (cf. 
Fig. [3]). However, below that point the data rise mono- 
tonically through 1 at yt ~ 3.1 and continue to rise for 
smaller yt, even below yt = 2 {pt = 0.5 GeV/c). The cen- 
trality trend near ?/( ~ 2 is closely (anti)correlated with 
the trend near yt ^ 5 (cf. Fig. [T^- left panel), strongly 
suggesting that the two widely-separated pt regions are 
physically connected by the parton energy-loss and frag- 
mentation process. More detailed centrality information 
for peripheral collisions is desirable to fill in the obvious 
gap (cf. Sec lVIHGI) . 

That the data for peripheral collisions are consistent 
with p-p data and unity confirms that references Snn 
and H]\[]sf for Au-Au collisions are properly defined. Spp 
was defined for p-p collisions as the limiting per-particle 
spectrum for rich [llj. By analogy, Snn should be 
consistent with the limiting Au-Au per-participant spec- 
trum for v 0. 

If a hydro contribution were present, we should expect 
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FIG. 10: Hard-component ratios for protons and five Au-Au 
centralities (thicker curves with changing line styles) relative 
to the N-N hard-component reference. The features are com- 
parable to Fig. |§1 

a deviation from the two-component model at the left 
side of the plot, shifting to the right with increasing cen- 
trality (v). The structure we observe appears to be part 
of the parton energy loss process. A real hydro contribu- 
tion should have a counterpart in the proton data, which 
we do not observe. 



B. Proton ratios 

Fig. 1101 shows Haa/ Hnn for proton data, with H^n 
defined in Sec. IVIBI The main differences from the pion 
data are the larger shift in the crossover point at unity 
relative to the model of Eq. and the large excess 
at yt ^ 3.5. This figure can be compared with v Raa 
trends in Fig. [7] (right panel). The excess in the proton 
hard component appears anomalous, but may be simply 
explained in terms of parton energy loss. 

The slopes of the proton data curves and correspond- 
ing model curves (dash-dot lines) match reasonably well 
as they pass through unity, although there is a horizontal 
shift of the crossing point (1.6 for protons, compared to 
0.4 for pions). The larger proton slope at unity (com- 
pared to pions) is expected in the hypothesis of Eq. (fT^]) 
because the width of the proton FD (and presumably 
the underlying proton FF for smaller parton energies) 
is ^ 1.5x smaller, and the rAA slope should therefore 
be ~ 2 times larger. The slope difference between pions 
and protons is thus quantitatively consistent with the FD 
width difference. 

The monotonic increase with decreasing yt begins near 



yt 4.5 for both pions and protons. The width difference 
between pion and proton FDs leads to quantitatively dif- 
ferent evolution at smaller yt (different slopes) but may 
refect the same energy-loss process for both species. The 
return of the proton ratio to unity at yt ~ 2.5 could 
have a counterpart in the pion spectrum at smaller yt, 
outside the pt acceptance. The return to unity below 
yt = 2.5, compared to the energy-loss model (sloped 
dash-dot lines), suggests that lower-energy partons lose 
relatively less (or no) energy compared to higher-energy 
partons. That possibility is discussed further in the next 
section. There is no apparent hydro structure in the pro- 
ton data. 



VIII. PARTON ENERGY LOSS AND COLOR 
SCREENING 

Raa data are modeled theoretically above pt = 6 
GeV/c {yt — 4.5) where, for central Au-Au collisions 
at 200 GeV, the constant value is near 0.2 for pions and 
protons. Within the restricted pt interval it is not clear 
whether that value derives from a reduction of the frag- 
ment yield, a leftward shift (negative boost) of the frag- 
ment distribution, a combination or some other process. 
Given the ambiguity I begin with a simple model of par- 
ton energy loss as a combination of limiting cases and 
consider possible deviations from that reference. 



A. Parton energy-loss reference 

A simple energy-loss model for changes in the hard 
component with A- A centrality is given by Eq. (jl3p , with 
attenuation {A < 1) and left shift (boost Ayt < 0) on 
yt ■ Approximating the logarithm of the hard-component 
ratio in Eq. (|14p with a Taylor expansion gives 

\og{HAA/HNN} ^ \og{A)~-Aytdlog{HNN)/dytil5) 

with ~d\og[HMN)/dyt = {yt-yt)/<yl^ or Uy^ (whichever 
is smaller) as obtained from the Gaussian-plus-tail model 
of FD Hmn in App. |^ Examples for A = 1 and Ayt val- 
ues corresponding to the five data centralities are plotted 
as the dash-dot line combinations in Figs. [91 and [TUl The 
sloped lines pass through unity at the Gaussian centroid 
yt = 2.66. The slopes are inversely proportional to the 
FD width squared. For the most central collisions in 
Fig. [5] alternative limiting case A ^ 0.21 with Ayt = 
is indicated by the labeled extension and seems not to 
describe the data. The negative-boost energy-loss model 
{A — 1, Ayt < 0) indicates the form of the ratio result- 
ing from shifting the entire FD down on rapidity with 
no other changes. The negative-boost model does not 
describe the data quantitatively, but its form provides a 
useful reference. 

The relative pt reduction corresponding to boost Ayt < 
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is given by 

p't « mo exp{2/t + Ayt}/2 « e^^* pt (16) 

From Figs. [51 and [TUl A-pt ~ —0.25 for pion and proton 
spectra above pt = 6 GeV/c in central Au-Au collisions, 
implying that about 25% of the parton (and fragment) 
momentum is "lost" in central Au-Au collisions for en- 
ergetic partons. It is not clear from the data whether 
parton energy is actually lost to a separate "medium" or 
instead is rearranged within the fragment distribution. 
The energy loss for lower-energy partons may be much 
less (cf. Sec. IVIII Pp . The seemingly large fragment at- 
tenuation at larger pt is apparently the consequence of a 
large power-law exponent and a small rapidity shift. 

B. vaa modeling and centrality evolution 



for pions and 

log[rAA{yui^)] = -Ayt{iy)d\og{H'^!^)/dytX (18) 
{tanh([yf - 3.24]/0.37) + l}/2 

for protons. Ayt(z^) is the negative boost for path-length 
(centrality) v, and H'j^j^ is with centroid shifted 

from 2.66 to 3.1 for pions and 4.2 for protons. The rAA 
models assume for simplicity that all centrality depen- 
dence is contained in multiplicative factor Ayt, the values 
shown as points in Fig. [12] (left panel). The small differ- 
ence between pion and proton rAA amplitudes at yt ^ 5 
may be due to the different exponents Uy^ = 5.5, 5.0 for 
pions and protons respectively. 

The ad hoc tanh factor in Eq. (|18p is well-determined 
by the proton data. However, the additional term in 
Eq. (jl7p is uncertain below yt ^ 2 where significant ad- 
ditional fragment yield appears. That uncertainty has 
implications for pion integrals and (pt) estimates. 



In Fig. 1 1 1 1 hard-component ratios rAA are plotted vs yt 
for pions (left panel) and protons (right panel) from Au- 
Au collisions in five centrality classes. The dash-dot lines 
represent the FD negative-boost energy-loss scenario — 
uniform shift of the entire hard-component FD to smaller 
rapidity — for central collisions. Despite the horizontal 
offset between data and model curves the similarities sug- 
gest that parton energy loss is dominated by a negative 
boost of some part of the underlying parton distribution. 
The offsets suggest that lower-energy partons experience 
less energy loss. 
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FIG. 11: Left panel: Hard-component ratios vaa for pions 
from Au-Au collisions for five centralities (solid curves). The 
dash-dot lines represent a rapidity shift Ayt (negative boost 
or red shift) of the entire N-N hard-component reference Hnn 
for 0-12% central collisions. The dashed curves are rAA model 
functions defined in the text. Right panel: Hard-component 
ratio rAA for protons from Au-Au collisions (solid curve). The 
features are similar to the left panel. 

The dashed curves in Fig. [TT] which closely model the 
rAA data are defined by 



log[rAA{yt;v)] 



-Ayt{v)d\og{H'^^)ldyt 
23exp{-1.9 2/t} 



(17) 



C. Comparison with energy-loss theory 

Parton energy loss in RHIC collisions is of major the- 
oretical interest [H, [H, |23| • I select one example [13] to 
provide a comparison with this analysis. 




FIG. 12: Left panel: Negative boosts Aj/t for yt spectra from 
five centralities of Au-Au collisions at 200 GeV (points) and 
theoretical predictions of relative energy loss AE/E for two 
parton energies (curves). Right panel: Same-side minijet peak 
widths from four centralities of Au-Au collisions at 130 GeV 
compared to values for p-p collisions {v — 1, hatched regions). 

In Fig. [T2| (left panel) I show values of Ayt obtained 
from this analysis vs participant path length ly. As noted 
in Sec. IVIII Ayt ~ Apt/pt ~ AE/E (signed numbers) 
for parton transverse momentum and energy. I also plot 
predicted relative energy loss AE/E for 10 GeV (solid 
curve) and 20 GeV (dashed curve) partons vs centrality. 

Two issues arise from this comparison. First, the 
data trend appears to exhibit a step-wise transition near 
ly = 2.5, already apparent in Figs. [9| and [TOl The 
data above the transition vary approximately as — Ayt ~ 
0.05 {i/ + 0.5). The theoretical centrality dependence is 
« ^ipart [13 ■ Because v ^ {ripart/^Y^^ the prediction 
is therefore AE/E oc , roughly gluon density v) x 
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outgoing parton pathlength {^-^ ly) . The predicted magni- 
tudes are close to the data, but there is no discontinuity 
in the theory. A step structure is also seen in minijet 
correlation peak systematics. Fig. [T^ fright panel) shows 
same-side minijet peak widths from 130 GcV Au-Au col- 
lisions (pt E [0.15,2] GeV/c), in which the rj width in- 
creases dramatically from p-p values (hatched regions) 
above v = 3 

Second, the energy dependence of relative energy loss is 
described as monotonically decreasing above Ejet 10 
GeV (Fig. 2 of That trend would be indicated 

by return of taa to 1 with increasing yt, which is not 
observed within the data pt acceptance. At smaller yt 
the trend of r^A for both pions and protons is a return 
to (and through) unity, which could be interpreted as 
a decrease in relative parton energy loss, an increase in 
fragment number due to energy loss or both. The proton 
data especially seem to indicate that low-energy partons 
lose little or no energy (cf. Sec. |VITTD|. 

For reference, at the transition point v ^ 2.5 (average 
number of N-N collisions per participant pair) there are 
on average about 19 participant nucleon pairs and 48 bi- 
nary N-N collisions (given = 42 mb), compared to 
191 and 1136 respectively for 6 = Au-Au. The frac- 
tional centrality is 58% (100% = N-N) and & = 11 fm. 



D. Color screening and charge screening 

The structure of Figs. [5] and [TO] compared to the ref- 
erence suggest that partons with smaller energies experi- 
ence less fractional energy loss compared to larger ener- 
gies. For instance, the hard component of proton spectra 
near pt ~ 1 GeV/c is consistent with the p-p hard compo- 
nent, even for central Au-Au collisions, whereas the FD 
negative boost above 6 GeV/c is consistent with a 25% 
parton energy loss. That interpretation, if correct, sug- 
gests the presence of color screening for partons — partons 
become increasingly "white" at smaller energy scales. 

Electric charge screening is observed in the energy 
loss of charged particles, especially heavy ions passing 
through crystalline media in certain orientations (chan- 
neling). The mean ion charge depends on the ion speed 
in such a way that straggling is reduced: more energetic 
ions carry larger exposed charges and lose more energy 
per unit path length, and conversely. Charge attachment 
from the medium to slower ions is more probable, leading 
to projectile charge screening. The consequent reduction 
of energy straggling is described as "bunching" [2£| . A 
similar screening phenomenon in RHIC collisions could 
lead to the observed enhancement of minijet structure 
relative to the N-N reference. 

Energy-scale-dependent color screening is consistent 
with e^-e" jet multiplicity trends. The distinction be- 
tween gluon and quark jet multiplicities, nominally deter- 
mined by Casimir color factors Cp and Ca, vanishes be- 
low Q = 10 GeV (5 GeV partons) [l^. Even at - 90 
GeV the gluon-to-quark jet multiplicity ratio Ug/uq ^ 1.4 



is significantly below Ca/Cf = 2.25 [sO]- That value 
(QCD) and 1 (non-QCD) are used to estimate parton 
energy loss in [1^, where the effect on Raa of differ- 
ing color charge is predicted to increase with decreasing 
fragment momentum. 



IX. THE PROTON-TO-PION RATIO 

We can use the two-component spectrum model to 
study the proton-to-pion ratio, which has received con- 
siderable theoretical attention [1, i, Fig. M (left 
panel) summarizes the model functions for pions and pro- 
tons obtained from this analysis. The dotted curves are 
the fixed soft components Smn- The dash-dot curves are 
the hard components Hmn ior v = 1 (N-N collisions), 
and the solid curves are Haa for i/ ~ 6 (6 = Au-Au). 
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FIG. 13: Left panel: Model functions inferred from this analy- 
sis. Dotted curves are soft components for protons and pions. 
Dash-dot curves are hard-component N-N (v = 1) references. 
The solid curves are hard components for u = & generated 
from the hard-component N-N references plus parton energy- 
loss parameterizations rAA(yt',v)- Right panel: Proton-to- 
pion two-component full-spectrum ratios on pt- The plot can 
be compared to Fig. 4 of [13 ]. 

Using the expressions for rAA {yt i^) obtained in the 
previous section densities 2/npart PAAh — {Snn + 
i^TAAiuu i^)Hnn} are defined for the two hadron species. 
Spectrum ratios pAAproton/pAApion are plotted in the 
right panel, labeled "protons" for v = 1, 6 (solid curves). 
The ratios for antiprotons (dashed curves) are obtained 
(according to observed data properties) by changing Uy^ 
from 5.0 for protons to 5.5 (the pion value) and multi- 
plying the resulting proton density by 0.8 to obtain the 
antiproton density. Those results can be compared with 
Fig. 4 of [l^ . Correspondence with the measured ratios 
is very good. The model ratios from this analysis are 
slightly larger due to the spectrum renormalization. The 
proton-to-pion puzzle is thus apparently transformed to 
details of parton energy loss and modified fragmentation. 

The full-spectrum ratios share the property of Raa 
that they mix soft and hard spectrum components, sup- 
pressing details at smaller pt- The change of the spec- 
trum ratio with centrality is actually modest. There is at 
the peak (2 GeV/c) only a factor 2x increase for central 
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Au-Au collisions relative to N-N collisions. In the left 
panel the hard-component ratio for N-N collisions (dash- 
dot curves) is 1 already at yt — 2.66 [pt ^ 1 GeV/c), 
descends to 0.5 near 4 GeV/c and then rises through 
unity again due to the apparent difference in the under- 
lying parton spectra (that for protons being harder). For 
central Au-Au collisions we observe an excess of protons 
above the FD mode and an even larger excess of pions be- 
low the mode. The differential hard components provide 
a much more detailed story than full-spectrum ratios. 

Comparing to ratios of fragmentation functions (FFs) 
from e'^-e~ collisions introduces further confusion be- 
cause the spectrum hard components are fragment distri- 
butions (FDs) — integrals of FFs folded with parton spec- 
tra. Even given that distinction, theory fragmentation 
functions such as the KKP parameterization typically de- 
scribe only the 10% most energetic fragments in an FF 
and seriously diverge from fragmentation data below that 
point |l2| , where A- A spectrum issues are most complex. 
In general, fully-differential formats such as Figs. [H] and 
1101 provide a clearer comparison of hadron species. 



X. COALESCENCE AND RECOMBINATION 

Ratios TAA are an improvement on Raa because they 
isolate hard components from soft, but absolute magni- 
tudes are lost. Additional insight is gained by studying 
transport of absolute quantities, particularly transverse 
mass, to test measure conservation during fragmentation 
or coalescence (competing models of hadron production). 

The two-component reference is an initial-state model: 
some hadrons are produced by immediate fragmentation 
of participant nucleons, some from large-angle parton 
scattering out of those nucleons and subsequent frag- 
mentation. Coalescence models describe hadron produc- 
tion by combinations of pairs and triples of partons both 
within and between the two parton spectrum compo- 
nents. Additional production of particles, momenta or 
quantum numbers relative to the reference could indicate 
extraordinary transport from longitudinal to transverse 
phase space. 



A. Hard-component Differences 

In Fig. [14] particle number transport (left panels) and 
rrit transport (right panels) relative to two-component 
references are plotted. There are two significant devia- 
tions from the reference: 1) particle (and mt) suppression 
at larger and 2) particle (and mt) excess at smaller pt. 
Because suppression at larger pt involves a small number 
of particles it is not visually apparent for protons and is 
a subtle feature for pions in this format. The conven- 
tional explanation for 1) is parton energy loss; that for 
2) is coalescence/recombination. This analysis suggests 
that given their closely correlated centrality dependence 
the two features may have a common origin. 




FIG. 14: Left panels: Number transport of pions (upper 
panel) and protons (lower panel) for five Au-Au centralities 
relative to hard-component reference Hnn ■ Right panels: nit 
transport of pions (upper panel) and protons (lower panel) for 
five Au-Au centralities relative to hard-component reference 
Hnn- 

B. Coalescence and Recombination 

Recombination was invoked for RHIC data in response 
to anomalous values of the proton-to-pion ratio at in- 
termediate Pt& [2,6] GeV/c for central Au-Au colli- 
sions d, [^, [l3|. The theoretical approach accounts for 
reduced suppression of protons in that interval, compared 
to the region at larger pt described by pQCD, by parton 
coalescence or recombination. 

The recombination mechanism assumes the presence 
of two partonic components, a thermalized, flowing par- 
tonic medium (QGP) with approximately exponential 
parton spectrum dominating at smaller pt, and a spec- 
trum of hard-scattered partons described by a fragmen- 
tation function (e.g., KKP) or a power-law function dom- 
inating at larger pt- The magnitudes of the two parton 
spectrum components cross over within the intermediate 
Pt interval. A coupling of the two components is proposed 
through parton coalescence to form hadrons, resulting 
in additional hadron yield in the intermediate pt region 
compared to spectra from elementary collisions. Differ- 
ences in spectrum details lead to the p/tt ratio anomaly. 

C. Assumed spectrum shapes 

Given the present heavy-ion and previous p-p two- 
component analysis there are problems with the spec- 
trum components invoked in the coalescence model. As- 
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suming that the soft component in elementary coUisions 
(partons and hadrons) is ah'eady thermal (including ra- 
dial flow) is unjustified. The hard component measured 
in N-N collisions, a fragment distribution built on a 
true power-law function (cf. App. is different from 
the fragmentation function (e.g., KKP) or conventional 
power-law function (modified Levy distribution) assumed 
in the theoretical models. The differences at larger pt are 
very significant. Results based on assumed spectrum fea- 
tures may be qualitatively different from data. 

The p-p two-component spectrum decomposition [Tlj 
revealed two independent non-thermal components re- 
sulting from fragmentation in two orthogonal directions. 
The notion that the soft component is the result of ther- 
malization should be contrasted with the possibility that 
the transverse soft component in N-N collisions is hadron 
fragments from the longitudinal "hard" component — 
partons from the participant nucleon PDF. At what point 
on Au-Au collision centrality does that situation change? 

A critical issue for the ReCo model is the locations of 
the spectrum crossovers which determine where the en- 
hancement by coalescence should appear. In 8, 9, j3| the 
pion crossover is in 2-4 GeV/c and the proton crossover 
is in 4-6 GeV/c. The result for the p/tt ratio is a peak 
within 2-4 GeV/c as shown in Fig. [T3] (right panel). How- 
ever, in this analysis we find that contrary to the coales- 
cence picture the ratio trend results from a combination 
of different hard-component widths plus the proton ex- 
cess centered at 2 GeV/c (not 4-6 GeV/c). The general 
shape of the ratio distribution is already established for 
N-N collisions where thermalization and QGP are not at 
issue, and the ratio increase with centrality near 2 GeV/c 
is due entirely to the localized proton excess, apparently 
related to parton energy loss. 



D. Two-particle correlations 



partons not explained in current coalescence models. 

E. Summary 

To be competitive, recombination/coalescence models 
must describe all aspects of the differential hadron spec- 
tra, including both hard and soft spectrum components, 
not just integrals or selected ratios. Whereas there is 
an unusually large baryon/meson ratio in one pt inter- 
val, there is an unusually large meson/baryon ratio in 
another (not predicted by present ReCo models). Both 
must be described by theory. Also, such models should 
transition smoothly from in vacuo N-N collisions (where 
there is no "thermalization") to central Au-Au collisions 
(where there may or may not be any "thermalization" ) . 

The picture that emerges from this analysis is one 
in which parton fragmentation remains the dominant 
hadronization mechanism in all pt intervals, but inter- 
actions among partons become important for larger den- 
sities (centralities). The hard components are strongly 
modified, but still maintain their integrity as parton frag- 
ment distributions. The differences between pion and 
proton FDs are substantial, but may simply be due to 
quantitative differences in the underlying FFs. Recombi- 
nation and coalescence models should be tested against 
that detailed experimental context. 

XI. INTEGRAL SPECTRUM MEASURES 

The two-component model with VAAiyuv) from this 
analysis describes spectrum data to their statistical lim- 
its. The model functions can be integrated to determine 
the centrality systematics of individual components and 
total spectra and compared to minimum-bias multiplicity 
distributions. 



Two-particle correlations present additional challenges 
for recombination. Even for the soft component in 
peripheral Au-Au (N-N) collisions we observe strong 
angular correlations suggesting a one-to-two (parton- 
to-hadrons) process leading to strongly charge- and 
momentum-correlated hadron pairs. We also observe 
strongly-correlated hadrons in the hard component down 
to 0.35 GeV/c which exhibit all the characteristics of jet 
correlations arising from fragmentation and which join 
smoothly onto the systematics of trigger-particle jet cor- 
relations as the selected hadron pt is increased. 

Thermal partons are assumed to be uncorrelated in 
the recombination model Q , in which case it is not clear 
how to obtain the hadron correlations observed in the soft 
component of N-N collisions. For the hard component, 
hadrons are modeled as SS or SSS combinations of shower 
partons. Back-to-back hadron pairs in the away-side jet 
would require a correlated combination of four, five or 
six shower partons from two jets to generate the hadron 
correlations, requiring a degree of coordination among 



A. Spectrum models 

The two-component models of pion and proton spectra 
from this study (3D densities) are summarized by 

-^P- = ^^r^{5o. + 0.012j.rAA.(2/t;:/)i/o.}(19) 

Upart 1.012 

2 0.062 po ^^ 
Pp = 1 no {Jop + 0.118t/rAAp(yt;i^)-nop}, 

npart i.iiS 

with po — 2.5/27r. The unit-integral model functions and 
hard-component ratios taa are defined above. The taa 
represent all deviations from the two-component refer- 
ence. If we set rAA = 1 and integrate Eqs. ([TO]) on yt the 
two-component constant for pions is a: = 0.012, whereas 
the K-N description requires xkn = 0.1 for the hadron 
yield at 200 GeV [l3|. The apparent contradiction must 
be contained in the rAA , which is confirmed here quanti- 
tatively. 
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B. Spectrum integrals 

In Fig. [15] (left panel) I show the integral of the charged 
pion (both signs) spectrum [first line of Eq. ([TO))] evalu- 
ated at six values of v (middle solid curve and points). 
The symbols refer to the integrated densities. The dot- 
ted line is the integral of pion soft component j^^^ . The 
solid curve corresponds to the taa from data, including 
parton energy loss and resulting fragmentation changes. 
The dash-dot line 0.85 po[l + 0.012(i/ - 1)] is the inte- 
gral if TAA = 1 in the spectrum model, describing the 
extrapolation from N-N collisions. 

Parton energy loss is thus responsible for a factor 5x 
increase in pion fragments for central Au-Au collisions 
relative to the N-N extrapolation on v. As noted previ- 
ously, there is a substantial uncertainty in the pion tj^a 
below yt ^ 2 which produces a comparable systematic 
uncertainty in the hard-component integral for more cen- 
tral collisions. The relative uncertainty in the pion hard 
component may be as large as 50%. The corresponding 
relative error in the pion solid curve in Fig. [TSlis 20%. 
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FIG. 15: Left panel: Spectrum integrals for pion soft compo- 
nent (upper dotted line) and hard component J dyt yt Haatt 
which sum to the middle solid curve. The dash-dot line ex- 
trapolates the pion hard-component reference from N-N col- 
lisions. The total and soft-component kaon contributions 
(lower solid curve and dotted line) are rough estimates. The 
upper solid curve is the total yield tt + K + p. Right panel: 
Similar curves for integrated proton spectra. The dash-dot 
line extrapolates the proton hard-component reference from 
N-N collisions. The differences between solid and dash-dot 
curves for pions and protons result from parton energy-loss 
factors Taa 1- 

In Fig. [15] (right panel) I show equivalent results for 
the proton spectrum (particle plus antiparticle) model. 
The proton hard component in N-N collisions is a much 
larger fraction of the total proton yield (11% compared 
to 1.2% for pions). The increase of fragment yield from 
parton energy loss is only 30% for protons (from dash-dot 
to solid curve) compared to factor 5x for pions. 

The lowest curves in the left panel provide a rough 
estimate of the kaon yield and centrality dependence. 
Adding kaons to the pion and proton results gives the 
top solid curve in the left panel, which is in fair agree- 
ment with the K-N two-component expression for the 
per-participant density of unidentified hadrons at 200 



GeV, 2.5 [1 -I- 0.1 - 1)] (top-most dashed line). We 
can now understand the large hard-component excess in 
Au-Au collisions over the N-N extrapolation. The main 
effect is the substantial increase in pion yield at small pt 
resulting from parton energy loss. The parton fragment 
yield for central Au-Au collisions is apparently 1/3 of the 
total multiplicity, and those fragments remain strongly 
correlated. 



C. (pt) systematics 

We can also use the model spectra to determine (pt), 
the ensemble mean pt defined by 

/>OC 

Po(j>t{v)) = / dptPt p{pt;i^) (20) 

= Ps{Pt)sNN +^Ph{PtM)HAA, 

with the densities scaled together appropritiately for each 
hadron species abundance. In the two-component model 
the mean pt also has two components, as in the second 
line of Eq. ([^0)1 : a soft component independent of cen- 
trality and a hard component. 
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FIG. 16: Left panel: Pion {pt) centrality trends from two- 
component model functions. The dotted line is the fixed 
soft component. The dashed line adds the extrapolated N-N 
reference vHnn- The solid curve with points is the com- 
plete model with rAA{jJt,v)- The open symbols are mea- 
sured means from extrapolated pt spectra [31.j . The dash-dot 
line is the two-component trend from p-p analysis assuming 
nch/rich ^ V [ll|. Right panel: Same as the left panel, but 
for proton spectra. 

Fig. [16] shows {pt) values obtained with the model 
functions from this analysis. The dotted lines indi- 
cate the constant soft component {pt)sNN- The dashed 
lines give the two-component trend expected for no par- 
ton energy loss: taa ^ 1 and vHaa ^ vHjsin with 
constant {ptjHNN- The solid curves and points are 
from the full data model with energy loss and evolv- 
ing {pt{i^))HAA- Ps and ph are given in Eq. HU). The 
soft- and hard-component reference means are (ptjSNN — 
0.40 ± 0.02, 0.76 ± 0.03 GeV/c and (]5t)ff„iv = 1-19 ± 
0.01, 1.12 ± 0.01 GeV/c, for pions and protons respec- 
tively in each case. 
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For pions (left panel) the spectrum mean (solid curve) 
can drop below the dashed line (arrow) because the addi- 
tional pion yield for central collisions appears to the left 
of the Hnn mode, as shown in Fig. 21 However, the solid 
curve should be considered a lower limit (cf. Sec. IXII CI 
for further discussion). For protons (right panel) the solid 
curve rises above the dashed line (arrow) because the ad- 
ditional proton yield relative to the reference hard com- 
ponent appears to the right of the mode in Fig. [S] 

The dash-dot line (left panel) summarizes the p- 
p multiplicity dependence of {pt) from (llj assuming 
nch/fich J^- The results are consistent within the 0.02 
GeV/c systematic uncertainty in extrapolating the spec- 
trum soft component to zero common to all rich or v [ll| . 
Also shown are mean values extracted from STAR spec- 
tra (open symbols) extrapolated from more-limited pt in- 
tervals: [0.2,0.7] GeV/c for negative pions and [0.5,1.05] 
GeV/c for antiprotons [3lj . 

The apparent mass dependence of (pt) and increase 
with A-A centrality is commonly interpreted to result 
from radial flow in heavy ion collisions, described for in- 
stance by a blast-wave model. Flow velocities are then in- 
ferred from spectrum fits [sij . From the present analysis 
the major contribution to the mass dependence is the soft 
component which does not change with A-A centrality. 
Most of the mass dependence is already present in N-N 
collisions. Support for an inferred "radial-flow" velocity 
(/3) is derived from the nearly linear increase of minijet 
yields, peaked nearpf — 1 GeV/c, with v ~ (npart/2)^/^. 
That dependence leads to the characteristic {dn/dr]Y/^ 
dependence of some inferred {(3) centrality trends. 

XII. MODEL UNCERTAINTIES 

This spectrum decomposition is based on three as- 
sumptions. 1) Deviations from participant scaling are 
minimal at small pt or yt (the hard component is small 
there). That region can be used to normalize measured 
spectra. 2) The first-order deviation from participant 
scaling and soft component So{yt) is the reference hard 
component (x vH^^iyt). Thus, So{yt) and Ho{yt) model 
functions are determined by variations with v in the limit 
v —^ 0. 3) Any nonlinearity at larger v is represented by 
factor rAAiut] t^) in the hard component. Those assump- 
tions are expressed by Eq. ([5]), and the limit process is 
defined in principle by 

PhHoiyt) = Mm {2/npartdpAA{yui')/di^} (21) 
PsSoiyt) = Mm {2/npartPAA{yui^) - l^PhHo{yt)}■ 
lIl practice, the relation between spectra for two periph- 
eral centralities and/or p-p collisions is used for the ini- 
tial determination of 5*0 and Hq, and the two-component 
parameters are refined once the entire system with dif- 
ferential plots is established. For each hadron species in 
this analysis the model parameters are the five spectrum 



normalizations, two shape parameters for each of 5*0 and 
Hq plus exponential (power-law) parameter ny^ and the 
hard/soft relative abundance x — ph/ps- 



A. Spectrum normalization 

Complete a priori ignorance of measured spectrum 
normalization was assumed, based on conclusions in 
App. [B] Initial guesses were based on the K-N two- 
component description of particle production measured 
by minimum-bias multiplicity distributions. The 200 
GeV p-p spectrum is normalized to 2D density po = 
2.5/27r. Au-Au pt spectra are normalized to po for — I, 
and the relative normalizations for other u are initially 
determined by the K-N trend. Relative normalizations 
are refined after the two-component analysis is estab- 
lished. Consequences of normalization uncertainties are 
negligible, as argued below. 

While assumption 1) (the low pt/yt region is slow- 
varying with multiplicity or centrality) is true for hadrons 
in p-p and protons in Au-Au collisions, it is not true for 
pions in Au-Au collisions (because of parton energy loss). 
However, the pion soft component falls rapidly while the 
hard component rises rapidly in the region near yt = 2. 
A 1% change in spectrum normalization produces large 
changes in hard components below yt ~ 2.5 {Hnn is 
about 1% of Snn at yt ~ 2, cf. Fig. S]). Requiring that 
the inferred hard component vary slowly with yt in that 
region reduces the relative normalization uncertainty be- 
low 2% for yt > 2.5. 

The form of the proton spectrum on yt is such that 
normalization uncertainty plays a negligible role in model 
uncertainties. Below yt ^ 3 the spectra exhibit near-ideal 
two-component behavior. The soft component is nearly 
flat where the hard component has significant magnitude. 
Applying the K-N normalization to the proton spectra 
according to App. [B] already results in convergence to a 
single Snn for small yt- Asymptotic approach of spectra 
to the soft component as ^ 1.5 is clear in Fig. [SI 



B. Soft component 

In general, the procedure of Eq. (PT|) is appHed to de- 
fine SNNiyt) = Ps Soiyt) and HNN{yt) = PhHo{yt). For 
instance, spectra for ly = 1, 2 would form "simultane- 
ous equations" for the two functions. Functions Spph 
and Hpph were obtained previously for hadrons in p-p 
collisions where we observe ideal two-component behav- 
ior [ll| . Thus, for pions Snntt was initially assigned the 
parameters {As, T, n^) = (20.3,0.1445 GeV, 12.8) from 
Spph ■ The result was a small but significant disagreement 
between inferred hard component Hnn-k and Hpph- By 
reducing rig to 12.0 consistency was achieved. The re- 
quired change in was attributed to the species mix 
of the hadron spectrum from p-p collisions and lack of a 
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power-law tail in the hard-component model. The soft 
parameter uncertainties are then similar to those in [ll[. 

No initial values were assumed for proton Snnp or 
H^Np] model parameters were determined solely by the 
iterative process defined by Eq. (|2ip. The final soft pa- 
rameters from this analysis for pt G [0.5, 12] GeV/c are 
{As, T, Us) = (5.25,0.223 GeV, 17), consistent with the 
proton soft component inferred from spectrum data for 
Pt e [0.3,3] GeV/c [25] (cf. Sec. IX of [lH). 



C. Hard component 

The hard spectrum components and their ratios to ref- 
erence distributions are the main result of this analysis. 
Sources of relative error for the hard components are 1) 
normalization error relative to fixed soft reference Spf]\j, 
common to all centralities (substantial for pions, mainly 
below ut = 2.5), 2) errors in relative normalizations of 
different centralities (negligible in all cases), 3) errors 
in Snn shape parameters, mainly Us (changes the Haa 
shape above the mode). 

Referring to Fig. 01 reduction of the overall pion nor- 
malization by 1% produces a 50% reduction in the i> ^ 2 
(peripheral) centrality at t/t ^ 2.25. The change for 
v ^ 5.5 (central) is negligible. Normalization issues are 
most critical for peripheral centralities where the hard 
component is small. Associated changes at j/t 3 are 
neghgible. 

Changes in Us = 12 by ±0.5 can be compensated at 
larger yt by adjusting the normalization of the soft com- 
ponent by factor 1 =F 0.05. However, the structure at 
smaller yt is then strongly distorted. In particular, the 
agreement between i/ ^ 2 data and the reference (unity) 
in Fig. [9] is spoiled at small yt- 

Given an optimized soft component, increase of the 
Hnn centroid by 0.05 causes a 20% decrease of the v 
2 data above the centroid and 20% increase below the 
centroid in Fig.O A decrease in the width 0.445 0.435 
causes an 20% increase oi ly ^ 2 data above yt — 3- 

Systematic variations for protons are similar. The 60- 
80% (most peripheral) data provide a strong constraint 
on the fixed model parameters. Evolution with centrality 
in Figs. [9] and [To] is then a characteristic of the data, not 
of the fixed model. The system of curves in those figures 
can be shifted or distorted collectively to some degree 
by adjusting the soft model parameters, but not their 
relative spacings. 

The (pt) determination for pions is quite uncertain for 
more central collisions because the shape of vaa is not 
defined by data below yt ^ 2 or pt ^ 0.5 GeV/c, and 
the fragment yield in that region increases by a large fac- 
tor ('--^ 5 X ) with centrality. The effect of the additional 
fragment yield on the spectrum mean depends on its lo- 
cation relative to 0.4 GeV/c (the soft-component mean) 
which is not determined by data. Thus, in Fig. \W\ (left 
panel) the solid curve plus points should be interpreted 
as a lower-limit estimate, and the dashed line an upper- 



limit estimate. In contrast, the proton (pt) (right panel) 
relative to the soft mean is well-defined for all centralities. 



D. Further comments on uncertainties 

It can be argued that the two-component K-N normal- 
ization trends were imposed on the analysis at the begin- 
ning, so the results in Fig. [15] are circular. However, as 
discussed above the normalization constraints from the 
data themselves actually determine the final normaliza- 
tions. That the K-N trends are approximately in agree- 
ment is then a consistency check, an indication of small 
systematic error. The more important lesson of Fig. [15] 
is the demonstration of soft- and hard-component con- 
tributions to the spectrum integral, their variation with 
centrality, and the dramatic role of parton energy loss 
in particle (especially pion) production in central Au-Au 
collisions. 

Figs. [9] and [10] show the residual differences between 
data and reference which contain information about par- 
ton energy loss. The 20% systematic uncertainty esti- 
mates (upper limits) in the absolute values of the data 
ratios should be compared with excursions on yt in those 
figures from 1/5 to greater than 10, a 50-fold variation. 
The relative structure in those figures is therefore very 
significant compared to the model uncertainties. If the 
two-component model is viewed as a "fit" to the data by 
the reference then Figs. [9l and [TOl show all the fit residu- 
als and reveal details of parton energy loss for all pt / yt 
in the acceptance. 



XIII. DISCUSSION 

Conventional interpretation of pt spectra from central 
heavy ion collisions is predicated on the formation of a 
thermalized and flowing bulk medium which strongly dis- 
sipates energetic partons. The present study, a differen- 
tial analysis assuming only linear superposition of two 
functional forms, reveals a different picture with surpris- 
ing simplicity. The hard component, all of which may 
result from parton fragmentation, extends well below 1 
GeV/c and dominates the spectrum centrality evolution. 



A. Two-component analysis 

The supporting context for two-component analysis of 
Au-Au pt spectra is provided by soft and hard compo- 
nents in elementary collisions: e^-e~ fragmentation func- 
tions described over the entire fragment distribution , 
p-p two-component analysis ll[ and extensive correla- 
tion analysis revealing minijet structure [28l. [s^. [ssl. [s^ . 
The question for heavy ion collisions is how the two com- 
ponents evolve with A-A centrality, and what that evo- 
lution tells us about QCD. The present analysis reveals 
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new information on parton stopping and fragmentation 
and provides critical tests of current theoretical models. 



B. The soft component 

The fixed form of soft component 5*0 is formally de- 
fined in this and the p-p two-component analysis |ll| as 
the limiting spectrum shape for z/ ^ or a; ^ re- 
spectively. As such, there is no physical model invoked. 
Particle production arises from an unspecified soft pro- 
cess (longitudinal fragmentation of individual participant 
or "wounded" nucleons). By hypothesis, the soft compo- 
nent does not change form with A-A centrality. Does 
some of the structure in Figs. [5] and (TU] actually belong 
to a changing soft component? The structure at large yt 
is consistent with pQCD expectations, and the structure 
at smaller yt is strongly correlated with it for pions and 
protons. The results of data analysis appear consistent 
with two-component assumptions concerning 5*0. 



C. The hard component 

Hard-component model HNNiyt) has a simple fixed 
form consistent with pQCD for larger yt, with recent 
results from p-p spectra for smaller yt 11], and with 
e~^-e~ fragmentation functions 12]. It appears to be a 
minimum-bias fragment distribution, the folding of an 
unbiased parton spectrum with conditional fragmenta- 
tion functions, extending down to 0.3 GeV/c. 

The Gaussian parameters of the Au- Au pion hard com- 
ponent are little changed from those inferred from p- 
p hadron spectra. There was no a priori information 
about the proton hard component. The form was in- 
ferred in this analysis from the centrality (v) dependence 
of Au-Au spectra following the iterative method applied 
to rich dependence of p-p spectra in 11]. The Hnn dis- 
tributions are similar: Gaussians on yt with exponential 
tails, having the same centroids but different widths and 
QCD exponents. The proton width is about 2/3 the pion 
width. 



D. Parton energy loss and color screening 

Measurements of parton energy loss inferred from nu- 
clear modification of parton fragmentation in the form 
Raa are expected to reveal properties of the QCD 
medium produced in RHIC heavy ion collisions. Hard- 
component ratio rAA measured over the entire pt or yt 
acceptance provides much more information, given the 
limited range of validity of Raa- In contrast to a pic- 
ture of "jet quenching" as absorption of partons in an 
opaque medium, the present analysis suggests that no 
partons are "lost" in A-A collisions. Their manifesta- 
tion (in spectrum structure and correlations) is simply 



redistributed within the fragment momentum distribu- 
tion, and the fragment number increases. A high-pt trig- 
gered jet yield may be reduced by a factor of five within 
particular pt cuts, but additional fragments emerge else- 
where, still with jet- like correlation structure ^28, .32.]. 

Figs. [H] and [TO] provide the first clear indication from 
spectra of how parton energy is transported within the 
fragmentation process. Much of the difference between 
hard-component ratios rAA for pions and protons may 
derive from the differing widths of the elementary (N- 
N) conditional fragmentation functions. The common 
indication from both species is an apparent decrease of 
relative parton energy loss with smaller parton energy, 
suggesting screening of color charge. 



E. The Cronin effect 

The Cronin effect is a modification of the pt spectrum 
in p-A collisions relative to a Glauber linear superposition 
of p-p collisions fv?\ . The observed effect is measured by 
nuclear modification factor RpA and consists of an excess 
{RpA > 1) in the pt interval [2,6] GeV/c. The Cronin 
effect is conventionally modeled as initial-state parton 
multiple scattering leading to increased kt- However, it 
was suggested recently that the Cronin effect could be a 
final-state effect due to recombination [s^ . 

Cronin enhancement has been contrasted with sup- 
pression of ratio Raa observed at RHIC for pt > Q GeV/c 
and attributed to jet quenching. Enhancement or sup- 
pression is measured relative to reference value i? = 1. 
What we learn from Figs. [3] and [7] is that the true refer- 
ence for such ratios is not unity. What has been called 
suppression may well be enhancement relative to the cor- 
rect reference, and conversely. In Figs. [7] and [TO] we see 
a large proton enhancement for central Au-Au collisions 
in just the pt region {yt G [3, 4.5]) where the Cronin en- 
hancement is observed, and with a similar form. It is 
possible therefore that at least part of the Cronin effect 
is a final-state phenomenon related to parton energy loss, 
albeit in cold nuclear matter. 



F. yt spectra and minijets 

The conventional RHIC picture of thermalized Au-Au 
collisions at 200 GeV is contradicted by the abundance 
of surviving minijets (fragments from a minimum-bias 
parton spectrum) observed in two-particle angular and 
momentum correlations [1^, [H, 13^. The "minijet ex- 
cess" directly confiicts with claims of near-ideal hydro- 
dynamics and "perfect liquid." Two-particle momentum 
correlations on TOt X mt {yt € [1.3,3]) apparently revealed 
the lower half of the pion Haa moving down on trans- 
verse momentum with increasing Au-Au centrality (3^ 
(cf . the systematics of Fig. H]) . Two-particle angular num- 
ber correlations reveal strong "minijet deformation" in 
the same pt range [28j . The energy and centrality sys- 
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tematics of two-particle angular pt correlations reveal the 
source mechanism of (pt) fluctuations and provide new 
details of parton- medium interactions [H, [sj] . 

The present analysis appears to resolve the contradic- 
tion. Although large-momentum partons do suffer sub- 
stantial energy loss (~ 25%), small- momentum partons 
apparently do not. Single-particle pt spectra in this anal- 
ysis are consistent with minijet manifestations in two- 
particle correlations which imply that low-energy partons 
are not absorbed by the medium. They are modified to 
some degree, but much jet structure remains at smaller 
Pt or yt- Color screening would explain the large abun- 
dance of minijet correlations observed in central Au-Au 
collisions. 



G. Further implications from this analysis 

Conventional spectrum analysis isolates different in- 
tervals of the Pt spectrum. Raa emphasizes larger pt 
and visually suppresses the strong minijet contribution 
at smaller pt , which is then reinterpreted in terms of soft 
physics models (e.g., radial flow). Spectra are fitted with 
monolithic functions motivated by thermodynamic, sta- 
tistical and hydro models whose parameters may not be 
meaningful but which are interpreted to reveal a struc- 
tureless, thermalized, flowing bulk medium. 

In contrast, the two-component spectrum model de- 
scribes several physical mechanisms consistently over the 
entire pt acceptance. Instead of a thermalized medium 
we find that 2/3 of the hadrons in central Au-Au colli- 
sions are part of a soft component apparently unchanged 
from N-N collisions, and 1/3 are part of a hard com- 
ponent which appears to be strongly-correlated parton 
fragments consistent with QCD expectations. In that 
context the "p/tt puzzle" is also a manifestation of par- 
ton energy loss. 

Since 1/3 of the hadrons originate effectively from 
rapidly- moving sources (parton fragmentation), the sig- 
nificance of statistical-model spectrum measures (chemi- 
cal potentials, decoupling temperatures, (pt)s) attributed 
to an expanding bulk medium can be strongly questioned. 
Upon close examination of pion and proton spectra no 
identifiable radial flow is apparent. Inference of an equa- 
tion of state or phase transition also seems problematic. 
The best evidence for a QCD medium may come from 
high-statistics minijet correlations [28l. [ssl. [3^ . 



XIV. SUMMARY 

In this paper I develop two-component model functions 
to describe pt and yt spectra for identified pious and pro- 
tons extending to 12 GeV/c. Spectra for five centralities 
from Au-Au collisions at 200 GeV are accurately repre- 
sented to the statistical limits of the data by a hard-|-soft 
model. This analysis extends a similar study of pt spectra 



for unidentified hadrons from p-p collisions at 200 GeV 
which provided reference data for the present analysis. 

The soft component has the form of a Levy distribution 
on transverse mass mt as in p-p collisions and appears to 
be independent of centrality. The hard-component refer- 
ence is a Gaussian on transverse rapidity yt, with expo- 
nential tail corresponding to the expected QCD power- 
law trend 1 jp^^ and required to describe pt spectra be- 
yond 6 GeV/c. The Gaussian width for protons is smaller 
than for pions by factor 2/3, but the peak modes are the 
same. 

The centrality dependence of the hard-component ref- 
erence is cx = '^ntiin/ripart, reflecting binary collision 
scaling. Evolution of the data hard component with cen- 
trality relative to the reference is represented by ratio 
rAA which generalizes nuclear modification factor Raa- 
Whereas the latter is a ratio of total pt spectra, includ- 
ing soft components, the former is a ratio of isolated 
hard components (data/reference). As such, taa reveals 
without distortion the centrality evolution of the com- 
plete parton fragment distribution over all pt or yt in the 
acceptance — all "parton energy loss" in Au-Au collisions. 

The evolution of rAA with centrality — the main result 
of this paper — is simply described. Parton energy loss 
produces a shift or negative boost Ayj of part of the 
minimum-bias fragment distribution (Gaussian plus tail). 
Much of the difference in taa structure for protons and 
pions may result from the substantial width difference of 
the N-N hard components. rAA structure at smaller yt 
suggests that energy loss of low-energy partons is small, 
and color screening may be important at smaller energy 
scales. Two-component spectrum analysis thus provides 
access to new aspects of parton energy loss in A-A col- 
lisions over the complete parton spectrum and resolves 
several open questions from minijet correlation studies. 

This work was supported in part by the Office of Sci- 
ence of the U.S. DoE under grant DE-FG03-97ER41020 



APPENDIX A: HARD-COMPONENT MODEL 

Pt spectra are expected to go asymptotically to power- 
law trend 1/p" at larger pt, reflecting an underlying 
power-law parton spectrum. In [llj | it was demonstrated 
that the hard component in p-p pt spectra is well de- 
scribed by a Gaussian on yt for pt < 6 GeV/c. However, 
that description is inadequate at larger pt, and a power- 
law trend must be restored to the model. The need is 
already apparent in Fig. 10 (left panel) of ^] where the 
data for rich = 11.5 rise above the soft Levy plus hard 
Gaussian model function for pt > 4 GeV/c. In this ap- 
pendix I describe how an exponential tail can be added 
to the Gaussian on yt to restore the QCD power law to 
the hard-component model. 
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1. Algebraic description 

The hard-component model function is defined on 
transverse rapidity i/f . ff the QCD power law is p(pt) cx: 
Pj""'' and 

mt pt 



then 



piyt) 

dlog[p{yt)]/dyt 



Vt 



-Pipt 



(Af) 



Since the region relevant to the power-law trend is j/t ~ 5, 
and systematic uncertainties in the exponent are compa- 
rable to 0.2, I define Uy^ = — 2 as the relevant expo- 
nential constant on yt. Since rih ^ 7.5 I expect Uy^ 5.5 
for data. 
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FIG. 17: Left panel: The logarithm derivative for a Gaus- 
sian (dashed line) and Gaussian plus exponential tail (solid 
lines). Right panel: Gaussian (dashed curve) and Gaussian 
plus exponential (power- law) tail (solid curve). The points 
correspond to yt values for the data used in this analysis. 

In Fig. [17] the algebraic strategy is illustrated. The 
hard-component Gaussian model from [ll[ (dashed curve 
g in the right panel) is 

g{yt) = Ah exp{-[iyt - yt)/ay,f/2}, (A2) 

with logarithm derivative d\og{g) / dyt = —{yt — yt)/<^y^ 
shown by the dashed line in the left panel. To add an 
exponential tail to the Gaussian the logarithm derivative 
in the left panel must be limited from below by fixed value 
— riyj (solid line in the left panel). The running integral 
of the resulting function (solid curve in the right panel) is 
then exponentiated to obtain the desired Gaussian with 
exponential tail as hard-component model i?o(yt). 



2. Numerical algorithm 

To illustrate a discrete numerical analysis I assume n 
arbitrarily-spaced data points at yt.i shown by the dots 
in Fig. [T7| (the yt values from this analysis) . Define ab- 
breviations Ig 



log(.g), IHq = log(i?o), Ig' = dlg/dyt 



{yt-yt)lal, lH^ = dlHo/dyt 
Obtain Igt = log[g{yt,t)] 
(right panel, points on dashed curve) 

Obtain forward differences dx : dxi^i 
with dxi = 0, for x = yt, Ig 



Obtain logarithm derivative Ig'i = dlgi/dyt^i 
Obtain forward average Ig'ij^i = Q-g'i^i + Ig'i)/'^ 
Obtain forward average yt.i+i = (yt,i+i + 2;t,j)/2 
Apply condition to produce exponential tail 



IH' = lg[ 



''Vt 



if lg[ > ~r 

if h/'< 



(A3) 



'at 



(left panel, points on solid lines) 
Obtain running integral 



(A4) 



i=l 



(right panel, points on solid curve). The hard-component 
model function with exponential tail is then Ho i — 
Ah cxp(ZiJo,i)> with Ah determined by the unit normal 
condition. 

Extension of the hard-component model to a Gaussian 
with exponential tail (power law on pt) results in qualita- 
tively better descriptions of data beyond pt = & GeV/c. 
Description of p-p and peripheral Au-Au collisions is very 
good to 12 GeV/c. The QCD power-law exponent is ex- 
pected to be Uh ~ 7.5, or Uy^ ~ 5.5. There is a possible 
difference between pion and proton exponents — 7.5 vs 7.0 
respectively — at a two-sigma level of significance. 



APPENDIX B: SPECTRUM NORMALIZATION 

Fig. [m shows published data from [13] without renor- 
malization (solid curves) compared to the final two- 
component model from this analysis (dashed and dot- 
ted curves). Published spectra were used in the form 
+ n^-, Up and np. When integrated, the spectra 
should lead to total yields whose centrality systemat- 
ics agree with minimum-bias distributions, a basic con- 
straint on spectrum normalization. 




FIG. 18: Published spectra (solid curves) for pious and 
protons for five Au-Au centralities with original normaliza- 
tion [l3|. The arrows indicate where the spectra cross over. 
Dashed and dotted curves show the final two-component ref- 
erence system from this analysis. 

The curve crossings at yt ^ 2.5 — 3 (arrows) suggest 
that there is a normalization problem. As published, the 
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data indicate that per-participant proton and pion yields 
below yt ~ 2.5 {pt ~ 1 GeV/c) decrease with increasing 
Au-Au centralifry. The figure can be compare with Fig. 3 
(left panel) of ll!| where a similar trend arises because 
the spectra for all hch have been deliberately normalized 
to unit integral. The crossing feature in that case resulted 
from the interplay between soft and hard components 
with changing control parameter rich- 



As an initial approximation the spectra were renormal- 
ized with two-component factors [l + 0.08(z^— 1)] for pions 
and [1 + 0.16(i/— 1)] for protons. The x values are based 
on Phenix data [36|. After the two-component analysis 
was established the factors were adjusted to final num- 
bers 1.07, 1.15, 1.23, 1.31 1.35 for pions and 1.05, 1.20, 
1.42, 1.68, 1.80 for protons based on spectrum structure 
near ~ 2, as indicated in the relevant figures. 
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